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Abstract 

Deformation properties of weakly bound nuclei are discussed in the deformed 
single-particle model. It is demonstrated that in the limit of a very small 
binding energy the valence particles in specific orbitals, characterized by a very 
small projection of single-particle angular momentum onto the symmetry axis 
of a nucleus, can give rise to the halo structure which is completely decoupled 
from the rest of the system. The quadrupole deformation of the resulting halo 
is completely determined by the intrinsic structure of a weakly bound orbital, 
irrespective of the shape of the core. 
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I. INTRODUCTION 



The nature of exotic nuclei with extreme isospin values is one of the most exciting 
challenges today, both experimentally and theoretically. Thanks to developments in experi- 
mental technology we are in the process of exploring the very limits of nuclear existence, 
namely the regions of the periodic chart in the neighborhood of the particle drip lines. The 
systems of interest are characterized by extreme values of isospin corresponding to large 
proton or neutron excess. 

On the neutron-rich side, there appears a region of loosely bound few-body systems, 
neutron halos (see Refs. for reviews). In these nuclei the weak neutron binding implies 
large spatial dimensions and the existence of the halo (i.e., a dramatic excess of neutrons at 
large distances). Theoretically, the weak binding and corresponding closeness of the particle 
continuum, together with the need for the explicit treatment of few-body dynamics, makes 
the subject of halos both extremely interesting and difficult. 

Neutron halos and heavy weakly bound neutron-rich nuclei offer an opportunity to study 
the wealth of phenomena associated with the closeness of the particle threshold: particle 
emission (ionization to the continuum) and characteristic behavior of cross sections 
existence of soft collective modes and low- lying transition strength [|8|- [l^] , and dramatic 



changes in shell structure and various nuclear properties in the sub-threshold regime IT^-|TB 



In this study we address the notion of shape deformations in halo nuclei. The importance 
of non-spherical intrinsic shapes in halo nuclei has been stressed in some papers, especially in 
the context of a one-neutron halo ^^Be and the nucleus ^B, alleged to be a one-proton halo. 
The ground state of ^^Be is a 1/2+ state. The low neutron separation energy, 5'„=504 keV, 
allows for only one bound excited level (l/2~ at 320 keV). The halo character of ^^Be has 
been confirmed by studies of reaction cross sections and the importance of deformation 
can be inferred from the large quadrupole moment of its core ^°Be, |(5|=229mb |18|. The 
halo character of *B has been suggested in Ref. where the large measured quadrupole 
moment of its I'^=2^ ground state, |(5|=68mb, has been attributed to the weak binding of 
the fifth proton in the lp3/2 state. (The existence of proton halo in ^B is still heavily debated 
For instance, Nakada and Otsuka in a shell- model calculation, demonstrated that 
the large value of \Q\ in ^B could be understood without introducing the proton halo.) 

The role of deformation in lowering the excitation energy of the 1/2'*' intruder level in 
^^Be has been recognized long ago. For instance, Bouten et al. |^1| pointed out that the 
position of the abnormal-parity intruder orbitals in odd p-shell nuclei can be dramatically 
lowered by deformation, and they performed the projected Hartree-Fock calculations for the 
parity doublet in ^^Be. In another paper, based on the cranked Nilsson model, Ragnarsson 
et al. demonstrated that the parity doublet could be naturally understood in terms 
of the [220] 1/2 ® 2si/2) Nilsson orbital. In particular, they calculated a very large 

triaxial deformation for the positive-parity level, and a less-deformed prolate shape for the 
negative-parity state. Muta and Otsuka p3| , p4| studied the structure of ^^Be and *B with 
a deformed Woods-Saxon potential considering quadrupole deformation as a free parameter 
adjusted to the data. Tanihata et al. |2^ concluded, based on a spherical one-body potential, 
that the positions of experimental drip lines are consistent with the spherical picture; they 
emphasized the effect of increased binding of the low-£ shell model states near the threshold 
that can give rise to the level inversion. In none of these above papers, however, have both 
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effects, i.e., the loose binding and self-consistency been simultaneously considered. 

Figure |l] shows the interaction radii for a series of Be isotopes deduced from measured 
interaction cross sections The relatively large radius for ^^Be has been interpreted 

as a sign for halo structure of this nucleus. It is, however, quite interesting to note that 
calculated deformations (as obtained in Nilsson-Strutinsky calculations) of the Be isotopes 
are found to vary in such a way that the corresponding nuclear radii reproduce the data 
quite well. In this case, no effects from the halo structure of ^^Be have been considered since 
the calculations are based on the modified oscillator potential. Although these calculations 
are somewhat unrealistic, the result displayed in Fig. |T] clearly stresses the importance of 
simultaneously considering both deformation and halo effects as, e.g., in ^^Be. 

Another, more microscopic, approach is the work by Kanada-En'yo et al. based on 
antisymmetrized molecular dynamics with improved asymptotics. They obtained very large 
quadrupole deformation for the 1/2+ state in ^^Be, and a less deformed l/2~ state. A similar 
conclusion has been drawn in the recent self-consistent Skyrme-Hartree-Fock calculations 

Also, very recently, the notion of deformation in ^^Be has been pursued by several authors 



p9|^3[ within several variants of the weak coupling scheme. In these models, the odd neutron 
moving in a Woods-Saxon potential is weakly coupled to the deformed core of ^°Be, and 
interacts with the core through the quadrupole field. The deformation of ^^Be is not treated 
self-consistently; either the strength of the quadrupole coupling is adjusted to the data to 
reproduce the quadrupole moment of ^°Be, or the deformation is adjusted to reproduce the 
energies of the 7=1/2 doublet. The advantage of the weak coupling approach is that the 
total wave functions are eigenstates of the total angular momentum and they have correct 
asymptotic behavior. In this context, a nice molecular analogy, in which the adiabatic 
coupling of the valence particle to the deformed core is applied, is the weak binding of an 



electron to a rotationally excited dipolar system ||3J], or to a neutral symmetric molecule 
with a non-zero quadrupole moment | |35| ]. 

In this study, we address the question of deformed two-body halos by considering the 
single-particle motion in the axial spheroidal square well. The corresponding Schrodinger 
equation can be separated into three ordinary differential equations in the spheroidal co- 
ordinate system. The properties of the deformed single-particle states, especially in the 
subthreshold region, are analyzed by making the multipole decomposition in the spheri- 
cal partial waves with well-defined orbital angular momentum. Large spatial extensions of 
deformed halos are discussed in terms of spherical radial form-factors. 

The paper is organized as follows. Section |I| contains the discussion of the generic 
properties of deformed halos. The method of evaluation of the prolate spheroidal wave 



functions is described in Sec. Ill, and the results of calculations are discussed in Sec. IV 



Finally, Sec. contains the main conclusions of the paper. 



II. DEFORMED NEUTRON HALOS, GENERAL PROPERTIES 



This section contains some general arguments regarding the concept of shape deformation 
in two-body halo systems. The material contained in this section is the extension of the 
analysis of Ref. ||3^ carried out for the spherical case. (For the corresponding discussion of 



three-body halo asymptotics, see Ref. 



3 



Let us assume that the weakly bound neutron moves in a deformed average potential 
U{r) (usually well approximated by a sum of the central potential and the spin-orbit po- 
tential). The neutron wave function can be obtained by solving the deformed single-particle 
Schrodinger equation 



Mr) = 0, (1) 



where K,^=y ~2me^/h'^ and is the corresponding single-particle energy (ej,<0). 

In the following considerations we assume the axial reflection-symmetric central poten- 
tial. (The generalization to the triaxial and/or reflection-asymmetric case is straightfor- 
ward.) Since the asymptotic properties of radial matrix elements depend very weakly on 
intrinsic spin, the spin-orbit term is neglected. Thanks to the axial symmetry, the single- 
particle states are labeled by means of A - the projection of the single-particle angular mo- 
mentum onto the symmetry axis (2;-axis) and parity, vr. Since the Hamiltonian considered 
is invariant with respect to the time-reversal symmetry, we shall only consider non-negative 
values of A. 

The deformed wave function can be decomposed into spherical partial waves with the 
well-defined orbital angular momentum £: 

^t{r)=Y.R^Ur)yiK{r), (2) 
where, due to the assumed reflection symmetry, 

TT = (-1)^. (3) 

At large distances {r>R, where R defines the arbitrary but fixed distance at which the 
nuclear interacion becomes unimportant), U{r) vanishes and the radial functions Rei,{r) 
satisfy the free-particle radial Schrodinger equation 



^ 2 d i[l + 1) 

dr'^ r dr ^ 



RiAr) = 0. (4) 



(Of course, Ri\y{r)=Re^y{r-) for r > R.) The solution of Eq. (||) is Riu{r) = Bihf{iK^r) 
where h'l is the spherical Hankel function and Bi is a constant, given by 

B, = R,,{R)/hj{iK,R). (5) 



The spatial properties of the system can be characterized by radial moments, 
and multipole moments {ilj^\r'^Yno\il)^) . Both quantities require the evaluation of radial 
matrix elements 



n.+2 p* 



(Mi/|r"|f Az/) = / r"+'Rl^^{r)R,Ur)dr = hu'Ku + (6) 







where / represents the contribution from the inner {r<R) region and O is the outer region 
contribution {r>R). Thanks to parity conservation [Eq. (^], i' = £ mod 2. The inner 
integral is, by definition, finite. The outer integral can be written as 
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= / hi*{tx)hj,{tx)x''+^dx. (8) 

J RKv 

In the limit of a very weak binding {n^ 0) one can replace the value of hliii^yR) in Eq. 
by the asymptotic expression valid for small arguments. This gives: 

Now, following the arguments of Ref. ||36|, one can demonstrate that for small values of n^, 
OnWu behaves asymptotically as 

0„«v oc — ^ { + const. ) , (10) 

where xq S> Rup is a small constant. Consequently, the asymptotic behavior of OnWu in the 
limit of small e strongly depends on quantum numbers n, i, and £'. Namely, for 

n + 1 : 0„«v diverges as {-e^f+^'-''-^'>^\ (11) 
n = £ + 1 : OnWu diverges as - ^ ln(-e,,), (12) 
n < i + £' — 1 : OnWu remains finite. (13) 



A. The normalization integral 



The norm of the deformed state [Eq (H)], N^y, can be expressed through the zeroth radial 
moment 

{M^uf = {^M) = E(^^'^kVA^/) = E (^o«A. + Omu) . (14) 

Consequently, according to Eq. (|ll]), the norm is divergent only if the deformed state contains 
an admixture of the s-wave. This is possible only for orbitals with 7r=+ and A=0. In this 
case, the norm behaves asymptotically as a/Oooo^oc(— e,^)~^/'', and the probability to find the 
neutron in the outer region, 

-Pouter = J _^ (-'-^) 

approaches one for zero binding. 
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B. The rms radius 



The root-mean-square radius of a deformed orbital is given by 

(AHr^lA^) = ^^J^^ = (ArA.)-^E(^AHr^|Mz.). (16) 

As discussed in Ref. with decreasing binding energy the integral O2UU diverges as 
(— C;^)^^/^ for i=0 and as (— e,y)^^/^ for i=l [see Eq. ([TTD ]. Therefore, in the deformed 
system, the rms radius diverges only if the Nilsson orbital in question contains a component 
of an s or a p state. This leaves only three classes of states for which the spatial extension 
can be arbitrary large: 

TT = +, A = : (r^) diverges as (— ej^)^"*^, 
TT = -, A = or 1 : (r^) diverges as {-e^y^^"^. (17) 

In the following, these states are referred to as "halo states" or "halos". Of course, this 
does not mean that other Nilsson orbitals cannot form very extended structures when their 
binding becomes very small. However, it is only for the states (0) that the rms radius 
becomes infinite asymptotically. 



C. The quadrupole moment 

The average quadrupole moment of a deformed orbital is given by 

{Au\r'Y,o\Au) = ^'^•^j''2''lf^ = (ArA,)-2^(Mz/|rVA^)(M|F2o|^'A), (18) 



where the angular matrix element is 




{iA\Y,o\i'A) = i/^W^-^(fA20|M)(r020|£0). (19) 



According to Eq. ([TT]), in the e^, — >■ limit the quadrupole matrix element diverges if £ + 
i' < 3. Since the quadrupole moment of an s state vanishes, the only diverging matrix 
element among the 7r=+ states comes from an s^d coupling. At small binding energies, 
the corresponding integral O202U behaves as (— ej^)"^/^. For negative-parity orbitals, the only 
diverging matrix element is the diagonal one, O2UU, which also behaves as {—e^)~^^'^. 

However, because of different asymptotic properties of the normalization integrals, the 
low-e behavior of the single-particle quadrupole moment of a weakly bound orbital does 
depend on parity. With e^— i>0, the quadrupole moment (|18D of the 7i=+ halo approaches 
the finite limit. On the other hand, for the 7?=- halos the norm of the state remains finite 
and the quadrupole moment behaves as (— Cj^)"^/^. 

It is instructive to consider the quadrupole deformation (32 extracted from the ratio 
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By splitting (r^y2o) and (r^) into contributions from the core (c) and from the valence (v) 
nucleons, one obtains 

~ T (H)e + (r2)v ■ ^^^^ 

For positive-parity halos (7r=+, A=0), the numerator in Eq. ( pij ) is finite while the 
denominator diverges as (— ej,)~^. Hence (32 is asymptotically linear in e^, i.e., it vanishes in 
the limit of zero binding: 

/?2(vr = +,A = 0)'-^°0. (22) 



On the other hand, for negative-parity halos (vr=-, A=0 or 1), the ratio (|2l|) is solely 
determined by the the p-wave components in the valence state: 

&(. = -.A)--^»|(lAK„|lA) = {+°-63 ifA = (23) 

That is, the deformation of the halo is solely determined by the spatial structure of the va- 
lence state wave function, independently of the shape of the core. The deformed core merely 
establishes the quantization axis of the system, important for determining the projection A. 



D. Higher moments and multipole deformations 



The above discussion is easily generalized to the case of higher multipoles. For instance, 
for n=4, the hexadecapole moment (r'^Kto) behaves asymptotically in the same manner as 
the quadrupole moment, i.e., it approaches the finite limit for the 7r=+ halos and diverges 
as (— Cjy)"^/^ for the 7r=- halos. However, the corresponding deformation P^, proportional 
to {r^Y^o) / {r^) approaches zero, regardless of parity of the halo orbital. 



III. THE MODEL 



In our study, the deformed potential U{r) has been approximated by a prolate spheroidal 
finite square well potential. Spheroidal infinite square well was early used by Moszkowski 
38| to discuss the properties of single-particle deformed orbitals. Merchant and Rae ||39[ 



investigated the continuum spectrum (e>0) of the spheroidal finite square well potential to 
calculate the particle decay widths of deformed nuclei. Since the main focus of our work 
is the behavior of bound single-particle orbitals very close to the e=0 threshold, particular 
attention was paid to a precise numerical solution of the Schrodinger equation in the limit 
of very small binding energies and/or large deformations. 



A. Prolate Spheroidal Coordinates and Parametrization of Nuclear Shape 

Assuming the z-axis to be a symmetry axis of the nucleus, the coordinate transformation 
between the prolate spheroidal coordinates (^, r], 0) and the cartesian coordinates {x,y,z) 
reads [H0H4 
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X 

y 

z 



av'(e-l)(l-r/2)cos</>, 



av'(e-l)(l-^^)sin, 



(24) 

(25) 
(26) 



where a > 0, 1 < ^ < oo, — 1 < "Z] < 1, and < < 2tt. 
Surfaces of constant ^=^o represent confocal ellipses, 



+ 



(27) 



with foci at (0,0, ±a), minor axis R±=a^J^Q — 1, and major axis R\\=a^o. (Since the purpose 
of this study is to investigate the generic features of weakly bound states in a deformed 
potential, the analysis is limited to prolate shapes. However, the calculations can easily be 
extended to the oblate side through a simple coordinate transformation.) 

It is seen from Eq. (27) that the parameter defines the shape deformation of a system. 
Indeed, for ^ 1, the surface (|27|) becomes that of a sphere with the radius a^o, while the 
limit ,^0^1 corresponds to a segment. Following Ref. we introduce the deformation 
parameter 6: 



2/3 



^2 
?0 



1/3 



.a -I, 



The volume-conservation condition [the volume inside the surface (|2 
on 6] yields 

Rn 



(28) 



should not depend 



(29) 



where Rq is the corresponding spherical radius. 

To find the relation between 6 and other quadrupole deformation parameters, one can 
compare the macroscopic quadrupole moment of the surface (^) 



Q2 



5 + 1 



with those obtained using other shape parametrizations 
between 5 and the oscillator deformation 5osc is 



'1 + l'^osc' 



2/3 



1 



-5 



(30) 



For example, the relation 



(31) 



For small values of 6osc, Eq- (|3T| ) gives 6 = |5osc- However, at a superdeformed shape 
(-^=2), both deformations are very close: 5osc=0.6 while 5=2^/^-1=0.587. Figure shows 
the family of shapes representing different deformations 6. 
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B. Bound States in the Prolate Spheroidal Well 



The deformed spheroidal square well potential is given by 



Uo for e < ^0 
fore>^o, 



(32) 



where Uq is the depth of the potential well (f/o<0) and depends on S through Eq. p8|). 

Expressed in prolate spheroidal coordinates, the time- independent Schrodinger equation 
(0) can be written as 



d 



d 



d 



d 



0. 



(33) 



Following Ref. ||4^, this equation can be separated into three ordinary differential equations 
by assuming ipi^^,!], (p) = _R(^)S'(?7)$(0). The functions R, S, and $ are solutions of the 
ordinary differential equations 



d 

A. 

drj 



d^ 



X 



Al 



2^dSM{c,ri) 



drj 



or 



2^2 



A2 



XaI-C 7] 



A2 



1 



i?A«(c,O = 0, 
SAi{c,ri) = 0, 



-A'$a(0), 



(34) 
(35) 
(36) 



where Xai is the separation constant and 



Cint = a^j2m{e -Uo)/fi for ^ < 
icext = ^a^J —2me/h for ^ > ^o- 



(37) 



In the following, Rxi{c,C,) and SAi{c,ri) are referred to as the radial and angular spheroidal 
functions, respectively. For positive values of e, scattering solutions for the spheroidal square 
well were solved in Ref. ||3^. Here, we concentrate on bound states with e < 0. 

The angular solution Sai^c,!]) can be expressed in terms of a series of the associated 
Legendre functions of the first kind 



(38) 



where the prime over the summation sign indicates that k=0, 2, ... if (/ — A) is even, and 
k=l, 3, ... if (/ — A) is odd (parity conservation). 

The radial functions Rai{c,^) are expanded in terms of spherical Bessel functions of the 
first kind, fi^\c^) = jk{cCj, and spherical Bessel functions of the third kind fk'\cC) = 
hk{cC) = jk{cC) + 'i^k{cCj- The internal radial function R^^{c,^) (^ < ^o) and the external 
radial function i?^/(c, ^) (^ > ^o) can be written as 
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where p=l or 3. 

Finally, the deformed single-particle wave function is given by 



^ ,N _ J ^nexcAii?y(Cint,O'5A«(Cmt,??)$A(0) for ^ < 

El 5„^xcA«^A/(«Ccxt,O'5Ai(«Ccxt,^?)*A(0) for^>^o 



where Qnt and Cext are defined in Eq. (|37D, and Uexc is the excitation quantum number 
labeling orbitals having the same quantum numbers A and 7r=(— 1)'. 

By matching the internal and external wave functions at ^=^0? one finds the eigenenergies 
e and the amplitudes An^^^Ai and B^^^^ai- The details of the calculation are outlined in 
Appendix ^ The procedure used to calculate the separation constant Aa; and coefficients 
(i^^(c) is discussed in Appendix 0. 



IV. RESULTS 

This section illuminates the general properties of deformed Nilsson orbitals discussed in 
Sec. H using the spheroidal square well potential. 

The single-particle energies of the finite spheroidal well with f/o=-80MeV and i?o=4fm 
are shown in Fig. ^ as functions of deformation 6. At a spherical shape the orbitals are 
characterized by means of spherical quantum numbers (ni). The deformed orbitals are 
labeled by parity tt, angular momentum projection A, and the excitation quantum number 
nexc which specifies the energetic order of a single-particle orbital in a given (7rA)-block, 
counting from the bottom of the well (e.g., nexc=l corresponds to the lowest state, nexc=2 is 
the second state, and so on). In the following, the deformed orbitals are labeled as [nexcAvr]. 
For example, the A=l orbital originating from the spherical shell Id is referred to as [11+] 
(see Fig. I). 



A. Radial Properties of Deformed Orbitals 

The dependence of the single-particle rms radius on binding energy is illustrated in Fig. ^ 
(spherical shape), and Figs. |^ and || (superdeformed shape). (In calculations, the binding 
energy was varied by changing the well depth Uq.) 

The spherical case has been discussed in detail in Ref. ||3^; here it is shown for the 
mere reference only. In all cases the asymptotic conditions [Eq. (|T^] are met rather quickly. 
Indeed, in the considered range of binding energy the values of (r^) for the Is state shown 
in Fig. ^ and the [10+], [20+], and [30+] orbitals of Fig. |^ approach an asymptotic limit 
[(— e)^"*^ dependence], and similar holds for the Ip state and [10-], [20-], [11-], and [21-] 
orbitals (see Fig. ^ which behave as (— e)^^/^ at low binding energy. The remaining states 
do not exhibit any halo effect, as expected. 

Figure |^ displays the probability Pouter [Eq- (HID] to find the neutron in the classically 
forbidden region, ^ > as a function of e for three superdeformed states with different 
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values of A. At low values of binding energy, the i=0 component completely dominates the 
structure of the [20+] state and Pouter 1- 

The radial form factors -R^Ancxcl^) appearing in the multipole decomposition [Eq. (^] 
carry information about the spatial extension of the wave function. They can be obtained 
by the angular integration: 



i?£An... (r) = J ^t^MYelir) dr. (41) 
Since in our calculations the total wave function ip^^^^^{r) is normalized to unity, the integral 

P,(Anexc)=/ |PMn_(r)|Vc?r (42) 
Jo 

represents the probability to find the partial wave £ in the state [nexcAvr]. Of course, 

^P,(Anexc) = l. (43) 

e 

Figures ^ and ^ display the radial formfactors for several orbitals in a superdeformed 
well assuming the subthreshold binding energy of e=-5 keV. For the 7r=+, A=0 orbitals 
(Fig. §), the i=0 component dominates at this extremely low binding energy, in spite of 
a very large deformation. Indeed, according to the discussion in Sect. [11 A|, the value of 



P^=o(Onexc) approaches one at small binding energies. In other words, the 7r=+, A=0 halos 
behave at low values of e like s waves. It is interesting to note that both the [20+] and [30+] 
orbitals are dominated by the 2s component; the corresponding £=0 form factors have only 
one node. For the 7r=- halo orbitals with A=0 and 1 (Fig. |^), the p component does not 
dominate the wave function completely (Sect. |II A|) , but a significant excess of a p wave at 
large distances is clearly seen. The radial decomposition of other orbitals (A > 1), shown in 
Figs. H and ^, very weakly depends on binding energy; it reflects the usual multipole mixing 
due to the deformed potential. 

The results presented in Figs. H and || illustrate the fact that the multipole decomposition 
of the deformed level depends on both deformation and the binding energy. Figures (7r=+) 
and (71"=-) show contour maps of P£(A?T,exc) for the A=0 orbitals as functions of e and 
5. The structure of the [10+] level, originating from the spherical Is state, is completely 
dominated by the i=0 component, even at very large deformations. A rather interesting 
pattern is seen in the diagram for the [20+] orbital originating from the spherical Id state. 
The i=2 component dominates at low and medium deformations, and the corresponding 
probability Pf=2 slowly decreases with 6, at large deformations approaching the (constant) 
asymptotic limit. However, a similar effect, namely the decrease of the i=2 component, is 
seen when approaching the e=0 threshold. In the language of the perturbation theory [0, 
this rapid transition comes from the coupling to the low-energy i=0 continuum; the i=0 
form factor of the [20+] orbital shows, at low values of e, a one-nodal structure characteristic 
of the 2s state (see Fig. At low deformations, the amplitude of the s component in the 
[20+] state is proportional to 5. Remembering that the norm, Eq. (p!^, behaves as (— e)~^/^ 
for the £=0 state, one can conclude that at low values of 6 and e the contours of constant 
Pf=o correspond to the power law 5^ oc (— e)^/^. This result, seen in Fig. tells us that the 
s component takes over very quickly even if deformation 6 is small. (A similar conclusion 
for the 1/2"*" ground state of ^^Be has been reached in Ref. ||33|| .) 
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The partial-wave probabilities calculated for the negative parity states [10-] and [20- 



presented in Fig. [IT] do not show such a dramatic rearrangement around the threshold. 
Namely, the [10-] orbital retains its i=l structure in the whole deformation region consid- 
ered, and the structure of the [20-] state at large deformations can be viewed in terms of a 
mixture of p and / waves. In the latter case there is a clear tendency to increase the i=l 
contribution at low binding energies, but this effect is much weaker compared to the [20+] 
case discussed above. 



B. Quadrupole Moments and Deformations 

Deformation properties of single-particle orbitals, namely the intrinsic quadrupole mo- 
ments (r^y2o) and quadrupole deformations P2, Eq- (PUP, are displayed in Figs. ^ and ^ as 
functions of binding energy at 5=0.6 (superdeformed shape). It is seen that the asymptotic 
limits discussed in Sec. |I| are reached at low values of e in practically all cases. In particular, 
the values of P2 for the 7i=+, A=0 orbitals approach zero with e ^ [Eq. (p2D], those for 
the 7r=-, A=0 states approach the limit of 0.63, and the value of /?2 for the [11—] orbital is 
close to the value of -0.31 [see Eq. (p3D]. The only exception is the [21-] orbital which at 
deformation 6=0.6 contains only 31% of the i=l component (see Fig. |^); hence the positive 
contribution of the If state to (r^y2o) still dominates. 

To illustrate the interplay between the core deformation and that of the valence particle 
[Eq. (PI )], we display in Figs. ^ and |13|: (i) quadrupole deformation P2 of the valence orbital, 
(ii) quadrupole deformation of the core, and (iii) total quadrupole deformation of the system. 
Here we assume that the core consists of all single-particle orbitals lying energetically below 
the valence orbital, and that each state (including the valence one) is occupied by two 
particles. It is convenient to rewrite Eq. (^) in the form: 

r, _ /?2,c(r^)c + /?2,v(r^)v _ P2,v + Xp2,c .... 



where 



X - f^- (45) 



For the halo states, x and /32,tot — > P2,v The results shown in Figs. |12]and[T^ nicely 
illustrate this behavior. Namely, for very small binding energy the total deformation of 
the system coincides with that of valence, regardless of the core deformation. These results 
illustrate the deformation decoupling of the deformed halo from the rest of the system. A nice 
example of this decoupling has been discussed by Muta and Otsuka [5^, who demonstrated 



that the halo proton in occupying the 1^3/2, A=l weakly bound orbit produces the oblate 
density distribution which greatly reduces the large quadrupole moment of the prolate core. 

C. Deformation Softness of Halo Systems in the Mean-Field Calculations 

The effect of decoupling of the valence particles from the deformed core in the limit of 
a very weak binding suggests that in such cases the constrained Hartree-Fock (CHF) or 
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Nilsson-Strutinsky (NS) calculations would produce very shallow potential energy surfaces. 
Indeed, in the CHF theory the nuclear Hamiltonian H is minimized under the constraint 
that the multipole operator that defines the intrinsic shape has a fixed expectation value 
{Q)=q. The intrinsic wave functions are found by minimizing the Routhian 

H' = H- PQ, (46) 

where /3 is the corresponding Lagrange multiplier. If Q is the quadrupole moment and the 
nucleus is weakly bound, then, especially in the case of halo systems, (Q) is very sensitive to 
small variations in the single-particle energy e,^ of the last occupied single-particle orbital. 
In particular, for the 7r=- halos, (Q) can take practically any value without changing the 
HF energy (H). This means that the numerical procedure used for searching for the self- 
consistent solution can be rather susceptible to uncontrolled variations of q with e,^. 

In the Nilsson-Strutinsky calculations, the bulk part of the binding energy comes from the 
the macroscopic energy -Emacro- Commonly used is the Yukawa-plus-exponential macroscopic 
energy formula MM which accounts for the surface thickness. The corresponding generalized 



surface energy reads 

where is the surface-energy coefficient, RQ=rQA^^^ , a is the surface diffuseness parameter, 
a = \r — r'\, and V denotes the volume enclosed by the deformed nuclear surface. The 
latter has been defined in our study by means of the axial multipole expansion in terms of 
deformation parameters Px: 



R(n) = c{p)Ro 



1 + 



(48) 



with c(/9) being determined from the volume-conservation condition. 

As demonstrated in Ref. for small deformations the generalized surface energy, 
Eq. (^7|), is given to second order by 

E, = E,isph) + J2cx{Of3l (49) 

A 

where the expansion coefficients cx solely depend on the dimensionless parameter 

C = ^ = ^. (50) 



It can be shown that the function ca(C) becomes negative below the critical value of (c 



which is roughly proportional to the multipolarity A. Consequently, the generalized surface 
energy is stable to /3a if C > Cc ~ 0.8A, or 

aX < Ro/0.8. (51) 

According to Eq. (|5TD, for a given nucleus both large multipolarity and large diffuseness can 
trigger the shape instability. This conclusion also holds for the finite-range droplet model 
mass formula lEsl. 
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The weakly bound neutron-rich nuclei and halo systems are characterized by very diffused 
density distributions. For instance, it has been predicted in Ref. [l^ that the average 



diffuseness in neutron drip-line nuclei can increase by as much as 50% as compared to the 
standard value representative of nuclei around the beta stability line. The effect of the large 
diffuseness on the macroscopic energy is illustrated in Fig. O, which displays -Emacro with 



the parameters of Ref. |^ for a light ^4=20 nucleus as a function of deformations (32, /?4, 



and /Sq. Since at low deformations different multipolarities are decoupled, Eq. (^Uj), they can 
be varied separately and the remaining ones are put to zero. The calculations are performed 
for three values of a. It is seen that the general rule given by Eq. (^) holds. Namely, for 
larger values of a and A the macroscopic energy becomes unstable to shape deformation, 
mainly due to instability of Es (for very light nuclei the effect of the Coulomb term is much 
weaker). Interestingly, the effect is fairly pronounced even for quadrupole distortions; the 
potential energy curve becomes unstable to /?2 already for a=1.3astd- The above results 
indicate that in the microscopic-macroscopic approach both single-particle and macroscopic 
energy become extremely shallow to deformation for weakly bound systems. 



V. CONCLUSIONS 

In the limit of very weak binding, the geometric interpretation of shape deformation is 
lost. Consider, e.g., a deformed core with a prolate deformation and a weakly-bound halo 
neutron in a negative-parity orbital. According to the discussion above, the total quadrupole 
moment of the system diverges at the limit of vanishing binding (i.e., (r^Y2o) can take any 
value). On the other hand, depending on the geometry of the valence orbital, the total 
quadrupole deformation of the (core+valence) system is consistent with a superdeformed 
shape (vr=-, A=0 halo) or oblate shape (vr=-, A=l halo). For a 7i=+ halo, the quadrupole 
moment is finite but /?2 approaches zero. In the language of the self-consistent mean-field 
theory, this result reflects the extreme softness of the system to the quadrupole distortion. 



Figure 15 shows an example of such a situation: The two valence particles occupying the 
weakly bound [11-] orbital give rise to an oblate deformation of the system, in spite of 
the prolate deformation of the core and the prolate shape of the underlying spheroidal well 
(5=0.2). 

Shape deformation is an extremely powerful concept provided that the nuclear surface 
can be properly defined. However, for very diffused and spatially extended systems the 
geometric interpretation of multipole moments and deformations is lost. 

The presence of the spatially extended neutron halo gives rise to the presence of low- 
energy isovector modes in neutron-rich nuclei. The deformation decoupling of the halo 
implies that the nuclei close to the neutron drip line are excellent candidates for isovector 
quadrupole deformations, with different quadrupole deformations for protons and neutrons. 
Such nuclei are expected to have a very interesting rotational behavior and unusual magnetic 
properties. For instance, the rotational features of such systems (moments of inertia, B{E2) 
values, (7-factors) should be solely determined by the deformed core. 

An example of the above scenario has been predicted in the self-consistent calculations 
for the neutron-rich sulfur isotopes performed using Skyrme Hartree-Fock and relativistic 



mean field methods |47,48]. When approaching the neutron drip line, the calculated values 



of P2 for neutrons are systematically smaller than those of the proton distribution. This 
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example illustrates once again that in the drip-line nuclei, due to spatially extended wave 
functions, the "radial" contribution to the quadrupole moment might be as important as 
the "angular" part. 

Finally, it is interesting to note that the anisotropic (non-spherical) halo systems have 
been investigated in molecular physics. A direct molecular analogy of a quadrupole-deformed 
halo nucleus is the electron weakly bound by the quadrupole moment of the neutral sym- 
metric molecule such as CS2 PB[] . 
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APPENDIX A: THE COMPUTATION OF SINGLE-PARTICLE WAVE 

FUNCTIONS 



The amplitudes An^^^M and Bn^^^M in Eq. (^) can be found following the method outlined 
in Ref. ||39|. The matching conditions for the internal and external wave functions (^OD at 
,^ = ^0 lead to the following set of equations: 



I I 



(Al) 



and 



The matching conditions [Eqs. ( |AT| ) and ( [A2D ] should hold for any value of rj. To eliminate 
this degree of freedom, one can take advantage of the fact that the angular spheroidal 
functions Sm{c-, rj) form a complete orthogonal set in the interval — 1 < < 1. By multiplying 
both sides of Eqs. ( [Al| ) and ( [A2| ) by SXi'{iCe^t,v) ^-nd integrating over r], one obtains the 
matrix equation for An^^^M- 



J^An^^^AiMw = for alir. 



(A3) 



where 



Ma, = I 



-^Azk^Cext, ^0 



5o 



X (S'A/'(zCext)|5'A/(Cint))- 



(A4) 
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The eigenenergies en^^^K-w and the corresponding amphtudes An^^^M are found from the 
condition that det(M)=0. 



APPENDIX B: SEPARATION OF THE SCHRODINGER EQUATION 



In the hmit of very small binding energy and/or large deformations, particular attention 
should be paid to the numerical method employed for determining the separation constant 
Aa; [Eq. (p4|)] and amplitudes d^^ appearing in Eqs. (|38|) and (^9|). 

The separation constant Aa; and coefficients d^\c) can be obtained from the three-term 
recurrence formula [HOll: 



k{k-l] 



(2A; + 2A- l)(2A; + 2A-3) 
"2(A: + A)(A; + A + 1) -2A2- 1 
(2fc + 2A + 3)(2A; + 2A- 1) 



+ 



+ 



(A; + 2A + l)(fc + 2A + 2) 
(2A; + 2A + 3)(2A; + 2A + 5) 
{k + A){k + A + 1) - Xj,i 



.Ai 

"fe+2 



4' 



0. 



(Bl) 



It is interesting to note that for c — > [the last term in Eq. (pT|) dominates] l=£, the 
separation constant becomes Aa; = and the only nonvanishing coefficient corresponds 

to k=i-A. That is, the wave functions (^) and (|39| ) are those of the spherical well. Indeed, 
according to Eq. (^), the limit of c — >• is reached at the spherical shape (a — > 0). However, 
the limit of Cext — is also approached at deformed shapes if the binding energy is small. 
By defining 



ttfc = 4-a/ 



2A; + 1(A;- A)! 



2 (A; + A)!' 

the matrix equation ( |B1| ) can be reduced to the Hermitian eigenvalue problem ^9 



(B2) 



XAi{c)ak 




{k + A + l){k + A + 2){k- A + l){k- A + 2] 
(2A; + l)(2A; + 5) 



ttfc+2 



^2k{k + l) - 2A2 - 1 
(2fc + 3)(2A;- 1) 



2k -1\ 



{k + A- l){k + A)(A; - A - l){k - A) 
(2A;-3)(2A; + 1) 



(B3) 



Another way of computing Aa; and d^^ is to employ the recurrence relations for the ratio 

dk/ dk-2 



42,41 



d, 



k+2 

dk 



+ /3k + lk 



dk 



d. 



0, 



k-2, 



(B4) 



where (3k and 7fc can be expressed through factors appearing in Eq. (|BT|). 

The technique based on the recurrence relation ( |B4| ) is useful when calculating radial 
wave functions at large deformations. In general, the coefficients dk are largest in magni- 
tude at /c = / — A, and quickly decrease with k. On the other hand, the magnitudes of 
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spherical Bessel functions of the third kind, hk{ic^), rapidly increase with k. Conseqently, 
the product d^^cJh^+Ai'ic^) becomes numerically unstable at large values of k if dk and /ifc+A 
are computed separately. This problem can be overcome by writing the recurrence relation 
^41) in the form 



dk+2h 



k+A+2 



h 



k+A+2 



dkhk+A 



h 



k+A 



+ /3k + 7k 



dkhk+A 

dk-2hk+A~2 , 



and employing the property of Bessel functions hk{iz 

2k + 1 



hk+2 2k + 3 

hk z 



+ 



z 

2k - 1 



1 + 



hi 



h 



k+A 



h 



k+A-2, 



hk-2 . 



(B5) 



0. 



(B6) 
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FIGURES 



FIG. 1. Nuclear rms radii calculated in the modified harmonic oscillator model (squares) for 
a series of Be isotopes. They are compared to measured interaction radii (open circles). The 
calculated deviation from the smooth 1.1^^/^ behavior (dot-dash line) is due to deformation. 

FIG. 2. Shapes corresponding to quadrupole deformations S=0, 0.2, 0.6, and 1.2, assuming 
the volume conservation condition, Eq. (p9|). 



FIG. 3. Single-particle energies of the finite spheroidal well with UQ=-80M.eV and i?o=4fm 
as a function of deformation 5. At a spherical shape the orbitals are characterized by means of 
quantum numbers (ni). The deformed orbitals are labeled by parity (7r=+, solid line; 7t=-, dashed 
line) and by orbital angular momentum projection onto the symmetry axis (z-axis), A. 

FIG. 4. Dependence of the single-particle rms radius on binding energy e at 5=0 (spherical 
shape) and i?o=4fm. The potential depth is adjusted in each case to obtain the desired value of e. 
The values of (r^) for Is and Ip orbitals diverge at small binding energies according to the power 
law, Eq. (|17D. 

FIG. 5. Binding energy dependence of (r^) (solid line), (r^l2o) (dashed line), and P2 (dotted 
line) for several tt=+ orbitals at 5=0.6 (superdeformed shape). The asymptotic limits discussed in 
Sec. H are indicated. 

FIG. 6. Same as in Fig. ^ except for the tt=- states. 
FIG. 7. Probability Pouter, Eq. (15), to find the neutron in the classically forbidden region. 



^ > as a function of e for three superdeformed states [12+], [21-], and [20+]. 



FIG. 8. Radial form factors, Eq. (|4l|), for several 7r=+ orbitals in a superdeformed well (5=0.6, 
i?o=4fm) at the "subthreshold" binding energy e=-5 keV. The probability P^(Anexc)) Eq. (p2|) , is 
indicated. 



FIG. 9. Same as in Fig. except for 7r=- orbitals. 



FIG. 10. Contour maps of probabilities Pq and P2, Eq. (|2D, for the [10+] (top) and [20+] 
(bottom) levels as functions of deformation 5 and binding energy e. Calculations were performed 
for i?o=4fm. 

FIG. 11. Same as in Fig. |l^ except for probabilities Pi and P3 for the [10-] (top) and [20-] 
(bottom) levels. 
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FIG. 12. Quadrupole deformation /?2 of the valence orbital (dotted line), that of the core 
(dashed line), and the total quadrupole deformation of the system (solid line) as a function of the 
binding energy of the indicated valence orbital of positive parity. Calculations were performed for 
i?o=4fm and 6=0.6. 



FIG. 13. Same as in Fig. 12 except for 7r=- valence orbitals. 



FIG. 14. The macroscopic energy, i^macro, foi" Z=N=W (normalized to zero at spherical shape) 
given by the Yukawa-plus-exponential mass formula as a function of (32 (top), (3^ (middle), and 

(bottom) for three values of surface diffuseness: astd=0.68fm (solid line), a=1.3astd (dashed 
line), and a=1.5astd (dotted line). It is seen that the macroscopic energy develops deformation 
instability with both increasing diffuseness and the multipolarity of deformation. 

FIG. 15. Contour plot (in logarithmic scale) of the pr^ distribution in the (x, z)-plane. Here, p 
denotes the single-particle density, i.e., p{'r)=J2u-occ. The deformation of the spheroidal 

well is assumed to be 6=0.2 and the corresponding prolate shape is indicated by a dashed line. 
Left: distribution for a core consisting of iV=4 particles (counting from the bottom of the well). 
Right: pr^ distribution for the A^=6 system. The two valence orbitals occupy the weakly bound 
e=-5keV [11-] orbital. 
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